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PACS 03 . 65 . Ud - Entanglement and quantum nonlocality 

PACS 03.67.Mn - Entanglement production, characterization, and manipulation 
PACS 14 . 60 . Pq - Neutrino mass and mixing 

Abstract. - Flavor oscillations in elementary particle physics are related to multi-mode entan- 
glement of single-particle states. We show that mode entanglement can be expressed in terms of 
flavor transition probabilities, and therefore that single-particle entangled states acquire a precise 
operational characterization in the context of particle mixing. We treat in detail the physically 
relevant cases of two- and three-flavor neutrino oscillations, including the effective measure of CP 
violation. We discuss experimental schemes for the transfer of the quantum information encoded 
in single-neutrino states to spatially delocalized two-flavor charged lepton states, thus showing, at 
least in principle, that single-particle entangled states of neutrino mixing are legitimate physical 
resources for quantum information tasks. 
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Various branches of condensed matter, atomic physics, 
and quantum optics have evolved in recent years towards 
the investigation and development of schemes for quantum 
information and computation science [1] . To this aim, en- 
tanglement is a key ingredient and a crucial physical re- 
source. Different forms of entanglement have been proven 
to be equivalent to observable quantifiers of performance 
success in quantum information protocols either with dis- 
crete [2] or continuous variables [3]. In the present work, 
we investigate the operational meaning of entanglement 
in the context of elementary particles physics. We will 
show that single-particle (mode) entanglement associated 
to particle mixing can be expressed in terms of transition 
probabilities in flavor oscillations, and can be exploited for 
quantum informational tasks. 

The concept of single-particle entanglement has been 
introduced and elucidated in a recent series of important 
theoretical papers [4-6]. Its use has been discussed in 
various contexts of quantum information, including tele- 
portation, quantum cryptography, and violation of Bell 
inequalities [7,8], and later experimentally demonstrated 
with single-photon systems [9-11]. Existing schemes to 
probe non-locality in single-photon states have been sub- 
sequently generalized to include massive particles [12]. In 
the present work we extend the discussion to the arena 



of elementary particles and provide a general operational 
characterization of single-particle entanglement in this 
context by showing its connection with the transition 
probabilities in any quantum system oscillating between 
different modes. We then show how this form of entangle- 
ment is in principle a real physical resource for the realiza- 
tion of quantum information protocols by discussing ex- 
plicit experimental schemes for transferring it to spatially 
separated modes of stable leptonic particles. These results 
allow to place mode entanglement in neutrino oscillations 
on equal footing with mode entanglement in single-particle 
atomic and optical systems. Furthermore, we show how 
the single-particle entanglement quantifies CP violation 
in neutrino mixing. 

Flavor mixing of neutrinos for three generations is de- 
scribed by the 3x3 Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) unitary mixing matrix U(6,8) [13], where 
(9,6) = (9i2, 6»i3, 6*235 8) an( l Cij = cos%, Sij = sin%. 
The parameters Oij are the mixing angles, and S is the 
phase responsible for CP violation. Here, without loss 
of generality we consider only Dirac neutrinos. In the in- 
stance of Majorana neutrinos, two additional CP-violating 
phases are present, which, however, as it is well known, do 
not affect the physics of neutrino oscillations. The three- 
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flavor neutrino states are defined as 

\yW) = XJ(6,5) \v {m) ) (2) 

where W^) = (\v e ) ,\v ^ ,\v T )) T are the states with defi- 
nite flavor and \v [m) ) = (h), W2), W 3 )) T those with def- 
inite mass. Let us recall that both \v a ) (a = e,/i,r) and 
\i>j) (j = 1,2,3) are orthonormal, i.e. iv a \vp) = S a ^ and 
(Uj\v k ) = S jtk . 

Neutrino oscillations are due to neutrino mixing and 
neutrino mass differences. The neutrino states \uj) have 
definite masses rrij and definite energies Ej. Their prop- 
agation can be described by plane wave solutions of the 
form \vj{t)) = e~ lEjt \vj) . The time evolution of the flavor 
neutrino states is given by: 

\v U) {t)) = U(i)|£ (/) ), 

U(t) = \J(l5)U (t)\J{6,5)-\ (3) 

where W-* ) are the flavor states at t = 0, Uo(£) = 
diag{e- lE ^ ,e~ iE ^ ,e- lE3t ), and U(t = 0) = 1. At time t 
the transition probability for v a — > vp is 

Pu a ^ ff (t) = |(^|^W)| 2 = |U^(0| 2 , (4) 

where a, = e, fi,T . The transition probability P Ua ^ u „ (t) 
is a function of the energy differences AEjk = Ej — 
(j,k = 1,2,3) and of the mixing angles. Since the cur- 
rent neutrino experiments deal with ultra-relativistic neu- 
trinos, the standard adopted approximation is AEjk — 

A 2 

2 ^ k , where A.m^ k = — ml and E = is the energy 
of a massless neutrino (all massive neutrinos are assumed 
to have the same momentum jt). 

Flavor neutrinos arc identified via charged current 
weak interaction processes, together with their associated 
charged leptons. In the Standard Model (SM), where 
neutrinos are taken to be massless, flavor is strictly con- 
served in such processes. On the other hand, neutrino 
mixing, consisting in a mismatch between flavor and mass, 
is at the basis of neutrino oscillations and CP violation. 
The introduction of neutrino masses as a correction to the 
SM is a necessary condition to explain such effects. 

When neutrino mixing is considered, loop corrections 
produce violations of lepton flavor in the charged current 
vertices: however, these effects are extremely small and es- 
sentially vanish in the relativistic limit [14]. Consequently, 
neutrino states entering weak interaction processes, like 
the ones where flavor neutrinos are created or detected, 
must be eigenstates of flavor neutrino charges. The corre- 
sponding operators can be rigorously defined together with 
their eigenstates in the context of Quantum Field Theory 



(QFT) [15]. In the relativistic limit, the exact QFT flavor 
states reduce to the usual Pontecorvo flavor states, which 
define the flavor modes as legitimate and physically well- 
defined individual entities. Mode entanglement can thus 
be defined and studied in analogy with the static case [16]. 

Equipped with the above setting, one can establish 
the following correspondence with three-qubit states: 
We) = |l>* e |0)„JO>„ T> M = |0>„Jl)„jO> Vr , \y T ) = 
10)^10)^11)^. States \0) Va and \l) Ua correspond, respec- 
tively, to the absence and the presence of a neutrino in 
mode a. Entanglement is thus established among flavor 
modes, in a single-particle setting. Eq. (3) can then be 
recast as 

Wa{t)) = V a e(t)\l) v JO) v JO)^+V a ^t)\0) Ue \l) v JO) Ur 

+ U«rW|0},j0)„Jl) lv , (5) 

with the normalization condition |U a ^(t)| 2 = 1 

(a, (3 = e,/i, r). The time-evolved states W it)) are en ~ 
tangled superpositions of the three flavor eigenstates with 
time-dependent coefficients. It is important to remark 
that, although dealing with similar physical systems, both 
the framework and the aim of the present paper differ sub- 
stantially from those of Ref. [16]. In the latter, by exploit- 
ing the wave packet approach, the multipartite entangle- 
ment, associated with the multiqubit space of mass modes, 
has been analyzed in connection with the "dccohcrcncc" 
effects induced by free evolution. In the present work, 
by exploiting the plane- wave approximation, the entangle- 
ment is quantified with respect to the multiqubit space as- 
sociated with flavor modes, and is related to the quantum 
information encoded in the neutrino flavor states, which is 
in principle experimentally accessible, as we will show by 
devising an explicit scheme for the information transfer. 

States of the form Eq. (5) belong to the class of W 
states. These, together with the GHZ states, define the 
two possible sets of states with tripartite entanglement 
that arc incquivalent under local operations and classical 
communication [17]. In some instances, only two neutri- 
nos are significantly involved in the mixing. For example, 
only the transition <-> v r is relevant for atmospheric 
neutrinos, while only the transitions of the type v e <-> v a 
are relevant for solar neutrinos. For two- flavor mixing the 
mixing matrix U(#, 5) reduces to the 2x2 rotation matrix 
U(0), 

U(0)=( C ° S " Sin J V (6) 
w V — sin cos J w 

the evolution operator reads Uo(t) = diag(e~ tElt , e~ lE2t ) : 
and the time-evolved flavor states yield the Bell-like super- 
position (a = e, n): 

Wa{t)) = \J ae m)u e \0)u„ +V a p(t)\0) Ve \l) v „. (7) 
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Bipartite entanglement of pure states is unambiguously 
quantified by the von Neumann entropy or by any other 
monotonic function of the former [18]. Among entangle- 
ment monotones, the linear entropy has a special phys- 
ical significance because it is directly linked to the pu- 
rity of the reduced states, and enters in the fundamen- 
tal monogamy inequalities for distributed entanglement 
in the multipartite setting [18]. As one moves from the 
two- to the three-flavor scenario, multipartite entangle- 
ment measures are readily available in terms of functions 
of bipartite measures [19-21]. Representative of this type 
of measures is the global entanglement. It is defined as the 
sum of all the two-qubit entanglements between a single 
subsystem and each of the remaining ones [19], and can 
be expressed as the average subsystem linear entropy [20] . 
Global entanglement can then be generalized by construct- 
ing the set of mean linear entropies associated to all possi- 
ble bi-partitions of the entire system [21]. An alternative 
characterization of multipartite entanglement is given in 
Rcfs. [22,23]. 

Let p = IV'XV'I be the density operator corresponding 
to a pure state describing the system S partitioned 
into N parties. Consider the bipartition of the A-partite 
system S = {Si, S2, ■ ■ ■ , Sjy} in two subsystems Sa„ = 
{S' ll ,S' i2 ,...,S , .j n }, with 1 < i x < i 2 < ... < i„ < N 
(1 < n < N), and 5 Bjv _„ = {%,%,..., S jN _ n }, with 
1 < ji < h < ■ ■ ■ < 3N-n < N, and i q ^ j p . Let 

PA n = Ph,^,...,^ = Tr BN _ n [p] = Tr juh ,..., jN _ n [p] (8) 

denote the reduced density matrix of subsystem Sa„ after 
tracing over subsystem Sb n _„- The linear entropy associ- 
ated to such a bipartition is defined as 

S l A "-' Bn -\ P ) = ^T(l - Tr An [p 2 A J) , (9) 

where the d is the Hilbert-space dimension given by d = 
minjdim Sa„ , dim Sb n ^ 71 } = min{2™, 2 Ar- "}. Finally, we 
introduce the average linear entropy 

(st N - n) ( P ))= (^)" 1 E^" ;Bn - ) (p) ) do) 

where the sum is intended over all the possible bi- 
partitions of the system in two subsystems, respectively 
with n and N — n elements (1 < n < N) [21]. 

We can now compute the linear entropies (9) and (10) 
for the two-qubit Bell state \v a (t)) : i.e. Eq. (7), with den- 
sity matrix p( a > = (^(^(fa Wl- The linear entropy as- 
sociated to the reduced state after tracing over one mode 
(flavor) can be computed straightforwardly: 

S M =S M = 4|U Qe W| 2 |U QM (t)| 2 

= 4\U ae (t)\ 2 (l~\V ae (t)\ 2 ) 

= 4|U aM (t)| 2 (1 - |U QM (i)| 2 ) ■ (11) 

In Eq. (11) and in the following, we use the notation 
S L e ^ = S^'^ (p^), where the subscript a refers to the 



time-evolved state (channel), and the superscripts (e; p) 
refer to the considered modes (flavors). Clearly, for the 
two-flavor state (7), and in general for any two-qubit sys- 
tem, symmetry imposes = = (S^ 1 ^). Ex- 
pression (11) establishes that the linear entropy of the 
reduced state is equal to the product of the two-flavor 
transition probabilities. Moreover, for any reduced state 
p of a two-level system one has that Sl = 2[1 — Tr(p 2 )] = 
ADetp — 4Ai(l — Ai), where Ai is one of the two non- 
negative eigenvalues of p, and the relation Ai + A2 = 1 has 
been exploited. Comparing with Eq. (11), one sees that 
the transition probabilities coincide with the eigenvalues 
of the reduced state density matrix. 

In Fig. 1 we show the behavior of as a function 

of the scaled, dimcnsionlcss time T — . I . In the same 
figure, we also report the behavior of the transition proba- 
bilities P Uc ^ Vis and Pj, a _,„ . The plots have a clear physical 
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Fig. 1: (Color online) Linear entropy S L e ^ (full) as a function 
of the scaled time T = , 2g | . The mixine ansle 9 is fixed at the 

experimental value sin 2 6 = 0.314. The transition probabilities 
Pv e — tv e (dashed) and P„ e ->„ (dot-dashed) are reported as well 
for comparison. 

interpretation. At time T = 0, the entanglement is zero, 
the global state of the system is factorized, and the two 
flavors are not mixed. For T > 0, flavors start to oscil- 
late and the entanglement is maximal at largest mixing: 
P„ e ^„ e = P Ve ^, v = 0.5, and minimum at T = ir. 

We can now investigate three-flavor oscillations, and 
the associated three-qubit IF-like states (5). Similarly to 
Eq. (11), tracing, e. g., over mode r, one has 

4r ;T) = 4|U Qr (t)| 2 (|U Qe (t)| 2 + |U aM (i)| 2 ) 

= 4|U Qr (t)| 2 (l-|U QT (t)| 2 ). (12) 

The linear entropies for the two remaining bi-partitions 
are easily obtained by permuting the indexes e, p, r. The 
average linear entropy for the state (5) is then 

= ^|u Qe WI 2 |u Q ,WI 2 + |u Qe WI 2 |u Q .WI 2 

+ |U QM (i)| 2 |U QT (i)| 2 ). (13) 

In Fig. 2, we show S^ 13 '^ and (Sf^) as functions of the 
scaled time T = . The mixing angles By and the 
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squared mass differences are fixed at the most recent ex- 
perimental values reported in Ref. [24]. In order to track 
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is equal to the asymmetry under time reversal, defined as 



= P 



P 



(i) - Pu a ^uA-t) 



(14) 



In the three-flavor instance, such a quantity is different 
from zero for a nonvanishing phase 8. It is worth notic- 
ing that 53/3 ^cp = with a,/3 = e,fi,T. Introducing 
the "imbalances", i.e. the difference between the linear 
entropies and their time-reversed expressions: 



AS 



(a,/3;7) 
LA 



tf n \t)~s^\-t) 



(15) 



we can compute, e. g., AS^^ , 



and obtain: 



Fig. 2: (Color online) Linear entropies S i "'' 3 ' 7 ' and {S^ ) as 
functions of the scaled time T. Curves correspond to the par- 
tial linear entropies S^^ 1 '^ (long-dashed), S^ T (dashed), 



AS { L r ]T) = 4A^(|U eT (t)| 2 + |U Te (0| 2 - 1) 



(16) 



where the last factor is CP-even. In Fig. 4 we show the 
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(dot-dashed), and to the average linear entropy (S^' 1 ') 



behavior of the imbalances AS 



K/3;t) 



as functions of time, 



(full). Parameters Oij and Am% are fixed at the central exper- and see how they effectively measure CP violation, 
imental values [24]. 



the behavior of the entanglement, we plot in Fig. 3 the 
transition probabilities P„ c ^„ Q (a = e,/x, t). Comparing 
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Fig. 3: (Color online) Transition probabilities P„ e _n, a as func- 
tions of the scaled time T. Parameters Oij and Amfj are fixed 
at the central experimental values as in Fig. 2. Curves corre- 
spond to P Ve ^ Ue (long-dashed), P Vk -+ Vii (dashed), and P Ve ^>v r 
(dot-dashed). 



Fig. 2 and Fig. 3, we observe that, as one may expect, 
the more mixed are the flavors, the higher is the global 
multipartite entanglement of the system. Moreover, the 
partial linear entropies S^^ 1 ^ and S^ T '^ measuring the 
reduced bipartite entanglement, exhibit a similar behavior 
due to the strong correlation between the components 
and v T . As T > the probabilities P„ e _n, and P v& ^ Vt in- 
crease and oscillate while remaining close. Similar consid- 
erations hold for states 1^(0) an( i Wr(t)). Entanglement 
and flavor transition probabilities for these states exhibit 
very fast oscillating behaviors, related to the experimen- 
tally measured values of the mixing parameters. 



Because of CPT invariance, the CP asymmetry A 
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Fig. 4: (Color online) The imbalances AS' i "' /3 ' 7 ' as functions 
of the scaled time T. Curves correspond to AS^ 1 (long- 



dashed) and AS i e e ' r; ' l) (dot-dashed). The quantity AS^ T ' e > is 
vanishing. Parameters Oij and Arafj are fixed at the central 
experimental values as in Fig. 2. The CP-violating phase is 
set at the value 8 = n/2. 



In order to demonstrate that the form of single-particle 
entanglement encoded in the time-evolved flavor states 
\vy> (t)) is a real physical resource that can be legitimately 
used, at least in principle, for protocols of quantum infor- 
mation, we discuss an experimental scheme for the trans- 
fer of the flavor entanglement of a neutrino beam into that 
of a single-particle system with spatially separated modes. 
For simplicity, we will restrict the analysis to two flavors 
a = e, fx. Consider the elementary charged-current inter- 
action between a neutrino v a with flavor a and a nucleoli 
N [13]. The quasi-elastic scattering interaction yields the 
production of a lepton a~ and of an outgoing baryon X, 
according to the reaction: 



N 



X. 



(17) 



In the simplest instance, the nucleon N is a neutron and 
the baryon X is a proton p; the corresponding scheme is 
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Fig. 5: (Color online) Scheme for the generation of single- 
particle entangled lepton states. A neutrino beam is focused 
on a target where events of the following charged-current inter- 
action may occur: v a + n — > oT + p with a = e , fj,. The beam 
is assumed to have energy higher than the threshold value nec- 
essary for the creation of a muon. A spatially nonuniform mag- 
netic field B(r) is then applied to limit the momentum of the 
outgoing lepton within a certain solid angle Sli, and to ensure 
the spatial separation between the electron and muon spatial 
paths. The reaction produces a superposition of electronic and 
muonic spatially separated states. 

illustrated in Fig. 5. Given the initial Bell-like superpo- 
sition (7), the unitary process associated with the weak 
interaction (17) produces the superposition 

\a(t)) = A e |l) e |0) M + A M |0)e|l) M , (18) 

where A e | 2 + 1 A M | 2 = 1, and \k) a , with k = 0,1, represents 
the lepton qubit. The coefficients A Q are proportional 
to U Q| g(t) and to the cross sections associated with the 
creation of an electron or a muon. Comparing our single- 
lepton system with the single-photon system, the quantum 
uncertainty on "which path" of the photon at the output 
of an unbalanced beam splitter is replaced by the uncer- 
tainty on "which flavor''' of the produced lepton. The co- 
efficients A Q play the role of the transmissivity and of the 
reflectivity of the beam splitter. Moreover, by exploiting 
the mass difference between the two leptons, the desired 
spatial separation between the flavors can be achieved by 
applying a nonuniform magnetic field. It is also important 
to remark that the approach proposed in the present work 
can be applied even in extended neutrino models includ- 
ing one or more sterile neutrinos. In such cases, from a 
mathematical point of view the main difference is that one 
deals with more than three modes (flavors), while, from an 
operational point of view, the presence of sterile neutri- 
nos (undetectable to date) would introduce a mechanism 
of loss of quantum information by making the (observed) 
mixing matrix non unitary. Therefore we can conclude 
that, at least in principle, the quantum information en- 
coded in the neutrino flavor states can be transferred to a 
delocalized two-flavor lepton state, and the single-particle 
mode entanglement acquires an operational characteriza- 
tion that can be exploited for quantum information tasks 



using systems of elementary particle physics. 
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